Chern-Simons Inflation and Baryogenesis 
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We propose a model of inflation where the Chern-Simons interaction and vector fields play a central 
role in generating an inflationary epoch. As a result, in accord with the APS mechanism Q, the 
Sakharov conditions for baryogenesis are self-consistently satisfled, and we calculate the net baryon 
asymmetry index in terms of the time-dependent gauge configuration necessary for infiation, based 
on the chiral anomaly. Infiation begins with a large plasma density of interacting gauge fields and 
fermions, which interact through gravity and the Chern Simons term. The Chern-Simons term drives 
power from an initial white-noise spectrum of gauge fields into a narrow-band of superhorizon wave 
vectors. At the same time, the fermionic current and metric coupling amplifies the gauge field on 
superhorizon scales. This amplification of horizon sized gauge field produces the correct conditions 
to maintain more than 60 e-folds of inflation. Eventually the gauge fleld dissipates by producing 
the observed baryon asymmetry — ~ 10"^", through the chiral anomaly and inflation ends. 



I. INTRODUCTION 

The most successful paradigm of the early universe is 
cosmic inflation. Despite its utility in model building, a 
number of conceptual and technical challenges associated 
with inflation driven by fundamental scalar fields have 
been discussed in the literature Jj. Inflationary mod- 
els driven from vector fields have been considered in the 
past beginning with the work of Ford [2] followed by other 
authors [3, 4]. Building on these past investigations, we 
present a model of inflation with new ingredients involv- 
ing the interaction of CP-asymmetric abelian gauge fields 
and fermionic currents, as opposed to solely gauge fields, 
to generate a realistic inflationary epoch. 

Moreover, this mechanism naturally satisfies the 
Sakharov conditions for baryogenesis during infiation by 
dynamically intraconverting the initial configuration of 
the CP-asymmteric gauge fields into the baryon asym- 
metry of the universe, building previous work of [5|. To 
be concrete, we explore a model of infiation that is gen- 
erated from multiple copies of Abelian gauge fields in- 
teracting with fermions in an FLRW background. The 
model differs from previous attempts of gauge field in- 
flation since it does not rely on any explicit gauge mass 
term nor non-trivial self interactions of the gauge fields 
to drive infiation. 

We now wish to spell out the basic idea of the mecha- 
nism before its mathematical instantiation. In this model 
of inflation we find that the negative pressure equation 
of state responsible for infiation is generated by a nearly 
constant interaction energy between the gauge field and 
fermion current Vint ~ ^/jV'7''^ (^s opposed to potentials 
that are purely functionals of gauge fields V{Afj,A^'-)). We 
assume that the universe comprises of initial configu- 
ration of abelian gauge fields A° and a coherent state 
of fermion currents Jj]. In an expanding space-time, 
gauge field and fermion currents typically dilute. How- 
ever, we find new solutions of the coupled field equa- 



tions demonstrating that the gauge fields are amplified 
due to a fermionic and metric source term. Eventhough 
the fermion current dilutes with cosmic expansion, the 
gauge field amplification will compensate and a constant 
interaction energy between the fermion and gauge field 
continues to persist within the horizon scale, which drives 
inflation. 

It is well known that vector models of inflation gener- 
ically lead to large anisotropics due to the spatial orien- 
tation of gauge fields. However, anisotropics can be nat- 
urally suppressed by employing N copies of gauge fields 
that are randomly distributed on the initial hypersurface 
at the beginning of infiation. As a result the anisotropics 
are suppressed by a factor of -4=. 

A crucial role that the Chern-Simons term plays in 
this model concerns one-loop quantum effects between 
gauge fields and fermions in the standard model {i.e. the 
Chiral Anomaly). Consequently, the Chern-Simons in- 
teraction is active during the inflationary epoch so as to 
naturally satisfy all three Sakharov conditions and can 
naturally end inflation by the production of an observed 
baryon asymmetry when the gauge field is converted into 
leptons through the Chiral-Anomaly|5|]. Assuming an in- 
stantaneous reheating, we calculate the net baryon num- 
ber to be - ~ 10"^^° in terms of the gauge fields that 
sourced infiation, making a further connection between 
gauge field driven infiation and baryogenesis; a relation 
between the observed baryon asymmetry and the initial 
density of gauge fields necessary for the onset to infiation. 



II. THE THEORY 

In this model, we assume that the early universe is dom- 
inated by N copies of an abelian gauge field interact- 
ing with a fermionic current. The gauge field A'Z has a 
U{1)^ symmetry with an index, A= 1..N. Generic string 
compactifications can easily accomodate multiple copies 



of gauge fields. In particular, one can obtain embed- 
dings by stacking N coindicent D-branes which leads to a 
SU(N) gauge theory. Symmetry breaking patterns exist 
where the gauge theory breaks to an abelian subgroup 
SU{N) — >■ C/(l)^~^; we refer the reader to the review 
by Langacker and references within for more details [6|. 
These extra f/(l)'s are anomalous and the Chern-Simons 
term is necessary for gauge anomaly cancellation due to 
the Green- Schwarz mechanism Q [S^l- The effective field 
theory that we consider for inflation is |33| : 



S = Sd + d Xy/—g 

Jma 
1 



M^R 1 



Stt 



-^ue^^'e + m^e^ 



where Sd is the covariant Dirac action: 



•d^a 






q V'7 e^j^^%l)A^ 



We will consider the dynamics of the fermions for self- 
consistency in Appendix A. A^ 4 is a 4D space-time man- 
ifold whose signature is diag(— 1, 1, 1, 1). The tensor 
F^^, —d[^A\ is the field strength tensor of A"^. We de- 
note the sum over the indices labeling the U{\) copies as 
Tr[...], while q stands for a dimensionless coupling con- 
stant. The axial current for spinors is J'^ ^'ip'^^e'^'y^ip. 
Here ip and ^ are Dirac spinors, 7^ with / = 0, ...3, and 75 
are Dirac matrices. Finally we refer to Mp as the Planck 
energy, M* is the mass scale identified with the UV cutoff 
scale of the effective field theory, and 9 is responsible for 
CP violation. 

For the purpose of efficiency of the presentation we 
will evaluate the dynamics of 9 field in Appendix A and 
show that its energy density is twenty orders of magni- 
tude smaller than the gauge-fermion interaction, making 
it insignificant for driving inflation. In this model, it is 
the Chern-Simons and fermion current interaction with 
the gauge field that conspires to yield the negative pres- 
sure equation of state in order to initiate and drive infla- 
tion. To establish this, we will first compute the energy 
momentum tensor of the Lagrangian of eq^ : 



C=Tt 



Jaj9f3S 



pjSpaP^ 



4Af* 



FafiF 



■aP ^ 



~^KJ^ 



Using the relation T^^'' = — ^^^HI^, we find that the 
energy-momentum tensor is: 



T,'f=-(5'^/: + Tr 



F Fa n°'PnPt^—2— P a F a ^f^ — A T^ 



=Tr 



-S^ApJi+FJ^FZ-l5iig"''g^''F^pFp,-A,J,^ . (2) 



We will now demonstrate that the purely isotropic part 
oi A-J' = TT[Af^J^^] in the energy-momentum tensor dom- 
inates over the anisotropic terms for large values of N. 
The other anisotropic terms in the energy-momentum 



tensor will be suppressed for the following reasons. First, 
the anisotropy due to the F^^iF^ term will be negligible 
because the corresponding electric and magnetic fields 
redshifts as a^"* during inflation according to the RHS of 
eq. (|T3)) . Also, the anisotropy due to Tr[^(yJ'^J term 
can be suppressed by two mechanisms. First, if there are 
N copies of randomly oriented gauge fields and currents, 
then the correction to the energy momentum tensor will 
be 0,13: 



T,,. 



N 
A=l 



H'i^j) 



J A ■ J5] + 0(1)Va^A(, J,-) . (3) 



The second term arises from the stochastic random distri- 
bution in directions of the fields that are non- vanishing 
if « ^ j. Notice that the ratio between the anisotropic 
and isotropic part of the averaged stress tensor scales as 

® 1^ ~ ^A-j ^ Vn- ^°'" ""^^"^^^ of iV - 10^ we 
can treat the anisotropic part of the energy-momentum 
tensor as a perturbation. 

Finally, another isotropizing mechanism is to have a 
triad of gauge fields that are mutually orthogonal A^ — 
6^A{ri,x) [1]. For more details about orthogonal gauge 
fields we refer the reader to [8| . Note that we can incorpo- 
rate the orthogonality condition in our model if iV = 3. 
However, this particular value of N introduces a fine- 
tuning problem. Whereas SU{N) gauge theories in the 
large N limit are quite natural in string theory [9| . 

We now focus our attention on the isotropic part of 
the energy momentum tensor ([2]) which is g = {E"^ -[- 
i3^)/(2a^) +A-J. Where we have introduced the defini- 
tions of electric field E = A and magnetic field B = V x A. 
It is now clear that the energy density has two compo- 
nents. The first component is the electromagnetic part 
that redshifts as a~^, and the second term involves a 
gauge-fermion-current interaction that we are going to 
show to be roughly constant and responsible for infla- 
tion. We shall now demonstrate that the coupled field 
equations indeed yield an inflationary epoch if certain 
initial conditions are satisfied. 



III. THE INFLATIONARY EPOCH AND 
INITIAL CONDITIONS 

In what follows, we derive solutions to the field equations 
for the gauge field coupled to both the metric and the 
fermionic current. The key to generating inflation is that 
the scale factor exhibits inflationary behavior when the 
gauge field and the fermionic current interaction remain 
nearly constant during inflation due to backreaction of 
the gauge field. 

First, we seek to find solutions of the gauge field in the 
FLRW background by varying the action with respect to 
A, 

5'^^fd"'Kp-e^P^'''d^^95l^F^,/{AM,) + ^gJ^=Q. (4) 



In the Coloumb gauge A^ = (0,^4^), equation ((4]) yields 



A - V^A + V X A e/K'h - a^J5 = ., 



(5) 



We seek to find self consistent solutions by working in 
conformal coordinates and assuming the expansion of 
the Universe to be given by a quasi-de Sitter phase, 
namely 0(77) = ao/[l — H{rj — 770)] and 770 the initial 
time. In ([5]) ' = d/drj, while V^ — d/dx^ and J^ 
stands for the spatial part of the fermionic axial cur- 
rent. Without loss of generality, using circular polar- 
ization vector fields and setting A^^ to be vanishing, we 
can write the field equations in terms of Fourier modes 
A{x, r])= J d^k X;h Mv, k)h eh(k)e''^'^, where h = ±1 de- 
notes the two possible helicities. Our field equation for 
the gauge field then simplifies to 

Airj, k)h + k^Airj, k)h = -hk A{rj, k)J/Ah + a^Jh , (6) 

where the spatial fermionic current is given by 

V2Jh = Jul a = {i) (71 + hi-12) l5'4')/a ■ (7) 

Here, the spatial current Jh = ip (71 + hij2) V' is shown 
to be a constant in Appendix A. We can then rewrite (jS]) 
as 

i(?7, k)h + k^A{r], k)h = -hk A{i^, k)h9/M, + a^J,, . (8) 

The general solution for the left-handed gauge field is 
then[3i]. : 

A{r,, k)^=A"_ cosh{l3kv)+A"_ sinh(/3fc77)+S[iJ, /?, J_ , 77] , 

(9) 
where A° and A^ are determined by the initial con- 
ditions, J± — flpJi and the growth factor is /3^ — 
k{9/M^, — k). We find that the particular solution is: 



m.PJ-.V]^ 2H4l~H,) +^nx) 






with 



^(x) = e^ / 



X f e *dt 



(10) 



We will have exponentia lly g rowing fields provided that 
(3k is real (or k < ^/M») [3a]- The important behavior 
for the gauge field is in the particular solution which gets 
amplified as 0(77) due to the factor J_/[2i/^(l — Hr])]. 

The guage field, A{r],k)+, of opposite helicity also 
grows as 0(77) i.e. 

A(77,fc)+=A'|cos(7fc?7)+i^sin(7fe77) + E[iJ,7,J+,77],(ll) 

in which 7I — k(9 / A'U + fc), and the particular solution 

'^[H,i,J+,y] = S is 



^ 7J+ /I . , ■, f cos tdt , . f sin tdt , , „, 

E=i^|--fsin(y)/^^ + cos(y)/ —— ,(12) 



2iJ3 \y 



with y — [\ — Hr])^/H . Notice that the particular solu- 
tion (IT^ grows as 0(7;) due to the factor J^/[27?^(l — 
Htf)] for all values of k. Moreover, ([T2| also exhibits 

a peak around a resonant momentum k — 2m~j which 
arises from the Chern-Simon's interaction. 

Now that we have the Fourier modes of the gauge field, 
^{Vt^)± "^6 can integrate of all momenta to obtain its 
spatial dependence: A{x,r]) = /d'^k A(77,k)± e''*''. In the 
next section we will perform this integral explicitly in 
order to show that the term ASf is nearly constant. 



IV. CONSISTENCY OF INFLATIONARY 
DYMAMICS 

We now turn our attention to the Einstein Equations 
and seek a consistent infiationary solution 36] . The Gqo 
component of the Einstein Equations gives the first Fried- 
mann equation: 



SttG 



3^ - 

a4 - a4 



{E+E-+B+B-)- 
+ &TiG{A+J^ + A^J+). 



(13) 



The coupled system can be solved if the interaction 
term A- J is nearly spatio-temporarily constant during 
inflation. The fact that J+ and J^ dilutes as 1/a might 
raise concern, but notice that for both the gauge field 
components, the dominant term grows proportional to 
0(77), as is evident from ^ and (fTTj) 37] and dominates 
the r.h.s of ([T3| . If A J is constant we can solve the Fried- 
mann equation (J13p to find an infiationary scale factor. 



a{v) =ao/[l- H{ri~'qo)]. 



(14) 



where ao is a normalization factor such that 0(770) = 
ao, and we have defined the Hubble parameter to be 
MpH ~ \/{A ■ J)no'^o- We can easily obtain the comov- 
ing scale factor a{t) using the map drj — dt/a, which is 
a{t) = aosxp Ht, which gives rise to exponential growth 
in the scale factor. With the choice of initial conditions 
such that a{r])/ao = (1 — Hrf)^^ = expHt, the time co- 
ordinate is conformally mapped to the bounded range 

77elO,il-i). 

Returning to the Fricdmann equation (J13p we can 
check the validity of our ansatz for spatial constancy of 
A- J hy using the solutions of the gauge fields ([9l fTTj) 
|38| . Taking the Fourier transform of ^(77, k)± yields the 
energy density oiA-J: 

I-Kwi^ - HV){J"- [aI cos(7fc77) + SJF, 7^, J+, 77] 



-J 



A°_ cosh{l3kv) + ^[H,(3kJ-,v]\ }e*'^"=A-^,(15) 



where we have normalized the integration over k by the 
factor 2/C. Since we are working with an effective field 
theory, we must perform momentum integrations up to 
a cut-off /C ~ M^. If we define the momentum scale 



k — 9/M^,, we can express the Fourier transform in eq 
P^ in a compact way and recognize that k contributes 
a spatial dependence of the form exp[±ifcz]. Notice that 
the resuh of integrating over k yields the rapidly decreas- 
ing functions of time, fi{r/). Therefore, the energy den- 
sity (fTS]) becomes: 

J^A"_hirj) + :^Mrj)]. (16) 



The second term on the RHS of (|16l) are exponen- 
tially time decaying terms in the quasi de-Sitter back- 
ground and can be disregarded, since they will decay 
after the onset of inflation. We immediately see that 
when ]Cz < 1 the first term in (J16l) becomes a constant, 

Q J-J^ sin{Kz) 



,J^J- 



m 



Kz 



2 ^2^ , and dominates the energy- 
momentum tensor. This condition corresponds to vector 
modes inside the Horizon, z < 4-. In this regime the 
FLRW equations become: 



M^H^ 






(17) 



Recalling that J± — AqJ^, the above equation consis- 
tently leads to inflation provided that J-/H ~ J+/H ~ 
10-^M2, A^_:^A\< IQ-^Mp and oq = 10^. Also, as we 
will analyze in the appendix, self consistency with the 
field equations requires O/M^, =k :^ H :^ M* . With the 
above values for the initial gauge fields, we find values of 
the fermionic current to be J+~J^c:i 10^^^ M^, namely 
J+ijio) ~ J-ijio) ~ 10~^''M^. Furthermore, in Appendix 
B we analyze the gauge invariant vector perturbations 
induced by fluctuations of gauge fleld 5A{r], k) and show 
that they decay as ^. 

We have demonstrated that the solution of the gauge 
field in a time dependent background self-consistently 
leads to an energy-density A- J that is constant in the 
Hubble radius at the beginning of inflation. The infla- 
tionary solution is self consistent for an initial density of 
fermion current and gauge field that are constant in a 
region inside the initial Hubble radius. 



V. INFLATIONARY LEPTOGENESIS 

Primordial nucleosynthesis and the recent determina- 
tion of the cosmological parameters from the cosmic mi- 
crowave background observations by the WMAP satellite 



require an excess of baryon to entropy density ratio 
in the universe to be 



13| 



S 



.5 ±0.4) X 10 



-10 



(18) 



where ng — ni, — n^ and s is entopy density of radiation. 
In this section we argue that this infiationary mechanism 



provides the possibility of generating the baryon asym- 
metry during the inflationary epoch. The baryogenesis 
mechanism, while compelling, will need additional ingre- 
dients, such as a rapid reheating mechanism before we 
can claim with more confldence that baryogenesis actu- 
ally occurs. In what follows we present encouraging clues 
that the Sakharov conditions occur quite naturally dur- 
ing inflation. We will explicitly calculate the net lepton 
number produced by the gauge fields during inflation, 
using the one-loop ABJ anomaly. 

Over forty years ago, Sakharov stated the three neces- 
sary conditions to generate a matter-antimatter asymme- 
try dynamically from symmetric initial conditions |14| : 

1. The vertices in the model violates baryon/lepton 
number. 

2. CP should be violated in the baryon/lepton number 
vertices. 

3. CP and baryon number violating interactions 
should be active when the universe is out of thermal 
equilibrium. 

However, the Standard Model of Particle Physics (SM) 
cannot account for the observed baryon asymmetry be- 
cause baryon number violating interactions in the SM 
are loop suppressed [H, [ig|. Also, the only source of CP 
violation in the hadronic sector is in the Dirac phase of 
the CKM mixing matrix, which is too small to generate 
the baryon asymmetry observed today 



17| 



Assuming 

that the scale of inflation is larger than TeV, the out- 
of-equilibrium condition for baryogenesis can be created 
at phase transitions or through decay of massive parti- 
cles. The most attractive choice for a phase transition 
is that associated with electroweak symmetry breaking. 
However, that phase transition is probably not strongly 
flrst-order [l5|,ll6|. 

Since the 1980's, it has been known that the weak 
interactions contain processes, mediated by sphalerons 
{SU{2) instantons), which interconvert baryons and lep- 
tons and are thermally activated at temperatures greater 
than 1 TeV. Thus, we can also create the baryon asym- 
metry by creating a net lepton number at high tempera- 
ture throu gh o ut-of-equilibrium and CP-asymmetric pro- 
cesses [13, lis]- These types of scenarios are known as 
leptogenesis. In what follows, we demonstrate that the 
Sakharov conditions are interrelated due to the CP vio- 
lating gauge field that is responsible for sourcing infiation 
itself. We then calculate the net lepton asymmetry pro- 
duced by the end of inflation. 



A. CP Violation 

It was shown by [5| that inflation provides a natural arena 
for satisfying the three Sakharov conditions for baryoge- 
nesis 39] . In [5| , the lepton asymmetry was generated by 
CP violating gravitational waves during inflation. In this 



model, however, it is the gauge field rather than gravita- 
tional waves that is responsible for the lepton asymmetry. 
Specifically, the interaction between the pseudoscalar and 
Chern-Simons term in ([S]) sources CP asymmetric gauge 
field configurations. The pseudoscalar coupling creates 
left and right asymmetry in the circular polarized gauge 
fields, A+,A-. The gauge fields that are generated no 
longer have definite transformations under CP: 



[CP]A+ ^ A^ 



(19) 



This condition can explicitly be checked by noticing that 
the left and right-handed gauge fields © and (jlip have 
unequal amplitudes. A parity transformation exchanges 
the handedness of a gauge field plane wave but does not 
transform the amplitude. Note that this form of CP vi- 
olation is not explicit but dynamical since the birefrin- 
gent amplitudes are sourced by the coherent evolution of 
the pseudoscalar, 0, during inflation, which arises from 
soluitions of 9 presented in Appendix A. 



B. Lepton Number Violation 

In this model, lepton number is violated at one-loop 
through the chiral triangle anomaly in the standard 
model ^]. Through the ABJ anomaly, the gauge held 
intraconverts itself into a net fermion number during the 
inflationary epoch via the following equation: 



V 7^ 



^ aj3^ ^u 



.aPtii/ 



/(32^2 



where the leptonic current is: 



(20) 



(21) 



and we have assumed that the leptons couple only to 
the hypercharge sector, U{1)y of the minimal standard 
model. Since we have extended the gauge sector to in- 
clude N copies of gauge held, we are assuming that these 
abelian factors couple solely to leptonic charge in the 
minimal standard model. In general we have the freedom 
to couple the gauge flelds to other states that are not 
charged under the weak-hypercharge assignment (such 
as a dark sector). This non- minimal charge assignment 
may connect our model to a mechanism of dark matter 
production at the end of inflation [40]. While this other 
possibility is to be pursued in an independent work, we 
wish to make a few remarks on non-standard charge as- 
signments. 

It is interesting that the many copies of U{1) neces- 
sary for our model to generate isotropy can be related 
to having extra Z' boson(s), if we make a non- minimal 
assignment of the fermions to extra C/(l)'s then they can 
be associated as dark matter, which was considered by 
Brahm and Hall |10|. We can explicitly implement the 
coupling of leptons to one of U{1)y as a the weak hy- 
percharge current in the standard model and the other 
L/(l)" copies as a set of Z bosons. In this coupling, 



the standard model neutral current interactions of the 
fermions will be in standard form: 



^WC ^gJ'^ernA^ + ffi ^^iZ° 



iM 



(22) 



where g and gi are the SU(2) and U{1)y gauge cou- 
plings; A^ and Z^ are the photon and massive Z-boson. 
On the other hand, the extra C/(l)'s can couple to the lep- 
ton current via the identiflcation with extra Z' bosons. 
Following Langacker and references therein [6|, we can 
modify the leptonic couplings to the the extended inter- 
actions of an SU{2) x U{1)y x [/(!)'" (with n > 1) gauge 
sector: 



C^fc = gJ' 



em-^^ 



n+1 



(23) 



where gi.Z^^ and J^i are respectively the gauge cou- 
pling, boson and coupling of the standard model. While, 
ga , Z^^ are the gauge couplings and bosons for the addi- 
tional C/(l)'s. In this non-standard coupling we will have 
to consider the production of leptons via the decay of 
the Z' bosons, which is an issue that requires a detailed 
analysis on its own right. 

We expect the non-standard couplings of the extra 
gauge flelds to alter our mechanism. We speculate that if 
the couplings are the same order as that of the standard 
model, then dark fermions will also be produced due to 
the triangle loop anomaly in the form of eq. ([20)1 . It is 
then plausible that the non-standard couplings to a dark 
sector could provide a potential mechanism for dark mat- 
ter production by the end of inflation via a kinetic mixing 
between the standard model U{1)y and the C/(l)' sectors. 
Such scenarios called Secluded WIMP Dark Matter have 
been discussed by Pospelov et. al. [ll|, and it will be 



interesting to connect those models to our scenario in fu- 
ture work; although this is a topic outside the scope of 
our present work [41[ |42| . 

We now show that net lepton number resulting from 
eq (PU]) is non-vanishing precisely because the gauge held 
conflguration that is responsible for inflation is also CP 
violating. This relationship between two Sakharov con- 
ditions is unique to our model. We can immediately cal- 
culate the amount of matter asymmetry produced from 
the gauge fields that initiated infiation through eq. (|20l) . 
Through the triangle interaction the gauge-field converts 
itself into a net lepton number accumulated throughout 
the inflationary epoch, and (I^Hl) becomes: 



H- 



drj 



k (i(77, k)+A{Tj, fc)+-i(?7, k)^A{f], fc)_) . 

(24) 

From ([9|) and (jlip . using the time dependent solutions 
of the gauge field A{ri, k)± and of currents necessary for 
infiation, we find n'^10~^^M^. 

Wc can use the Friedmann equation to find the entropy 
density, s = l.Sg^n^, where g^ is the effective number of 
massless degrees of freedom {e.g. 100 for the MSSM) and 



n^ = 1.285* {HMpY/"^ . Hence, we calculate a realistic 
value for the baryon asymmetry index n/ s ~ 10~^°. Note 
that the time evolution of the gauge field, A± is crucial to 
obtain the correct lepton- number, which shows the inter- 
esting relation between the dymamics of the gauge field 
during inflation and leptogenesis. We intend to pursue 
the real time dymamics of the lepton number production 
in a forthcoming paper. 



C. Out of Equilibrium 

Lepton number and CP violation occur simultaneously 
during inflation through the gauge-fermion vertex and 
triangle loop diagram respectively. However, the expo- 
nentially expanding background space-time renders the 
universe to be far out of equilibrium with the fermion 
production. Furthermore, because we are generating the 
Sakharov conditions during inflation, this leptogenesis 
mechanism depends on the details of a reheating mech- 
anism to conserve the net-leptons produced during in- 
flation. We assumed a spontaneous reheating [25i] . how- 
ever if reheating occurs slowly to a reheating temperature 
Tr < T, then the baryon asymmetry index, — will be di- 
luted by a factor Tr/T. Therefore, to firmly establish 
that we have a successful baryogenesis mechanism will 
require a detailed analysis of reheating in this model; a 
topic that we leave for future work. 



VI. END OF INFLATION 

In most scalar models of inflation, inflation ends when 
the inflaton no longer satisfles the slow-roll conditions. 
Our model is similar in that the end of inflation is trig- 
gered by the late time oscillation of the 9 field about its 
minimum. Recall, that inflation is driven by a nearly 
constant interaction energy between the gauge-fields and 
fermionic currents, which dominates the energy momen- 
tum tensor. The solution of the field equations for the 
gauge fields responsible for inflation (eq [5]) arose from 
a phase where the velocity of the fleld is slowly-rolling 
(i.e. 9 ^ const): 

A{7^,k)hWA{r^,k)h = -hkA{f^,k)h9/M,+a^Jh. (25) 

Eventually, the 9 field, which couples to the vector 
fleld, evolves to the bottom of the its potential, be- 
gins to oscillate, and the slow-roll condition is violated, 

'^- 2 \ V(9) j -^■ 

Eventhough the scalar potential does not source infla- 
tion, its eventual oscillatory coupling to the gauge fleld 
will end inflation. Hence, following the work of Kofman 
et al. and Traschen et. al [2a,[23, at late times the scalar 
fleld approaches sinusoidal behavior: 



resulting in a Mathieu equation which yields a resonant 
particle production: 

i(?7, k)h + k'^A{r], k)h + $o sin(mi) ^(77, k)h + a^A = . 

(27) 
Notice that we can 



where the amplitude is $0 



kMp 
37r Mtf 



6 = 



Mp 



s[ii{mt) 



(26) 



use (f26| because of the fast roll condition (e > 1), as the 
kinetic term 9 dominates over Chern-Simons potential. 

As a result of the oscillation of the 9 fleld about its 
minimum, the gauge field ceases to have an exponen- 
tial time dependent growth to sustain inflation. This 
happens because the Chern-Simons coupling of 9 to the 
gauge fleld mimics a Mathieu equation (l?7|) for the gauge 
fleld. This leads to resonant particle production similar 
to what happens in parametric resonance in preheating. 
However according eqs. ([27|) and ([24| . CP asymmetric 
parametric resonance and the gauge fleld will decay into 
fermion number via the Chiral anomaly. As a result, the 
oscillation of the 9 fleld triggers a sequence of decay of 
gauge fleld into fermions that ends inflation. Currently 
one of us and fellow collaborators are pursuing a pre- 
heating mechanism based on the coupled dynamics of 
the fermion, gauge fleld and oscillating 9 fleld numeri- 
cally [H]. 



VII. CONCLUSION AND DISCUSSION. 

In this work we demonstrate that it is possible to obtain 
an epoch of cosmic inflation from abelian gauge flelds 
sourced by the Chern-Simons current. This interaction 
leads to an exponential time dependence of the gauge 
fleld that compensates the diluting fermionic current to 
give a nearly constant energy density pAj ^ A-J. Power 
is dumped into a narrow band of modes, which drives 
inflation. Furthermore, we demonstrate that the CP vio- 
lating gauge flelds responsible for inflation can get inter- 
converted into a net lepton number at the end of infla- 
tion. We calculate that this number gives a realistic lep- 
ton asymmetry. However, we assumed that the entropy 
production occurs from spontaneous reheating. We did 
not pursue a reheating mechanism inherent to this model. 
We have argued that parametric particle production of 
the gauge particles due to coherent oscillations of the 9 
fleld triggers the end of inflation, which could yield a pre- 
heating phase, and we leave this numerical analysis for 
future work. 

We plan to pursue the genericity of the assumptions of 
our initial conditions by implementing our mechanism in 
a cosmological fermionic bounce mechanism, which may 
set up the initial condition for the gauge and fermion 
fleld. Some of us are currently pursuing this ques- 
tion of generating an initial fermion current in the con- 
text of a bouncing cosmology where a fermionic current 
could bounce from a contracting FLRW to an expanding 
FLRW, by using techniques in Loop Quantum Cosmol- 
ogy developed in [29| . 

Finally, an explicit analysis of the generation of scale 



invariant spectrum of density perturbation is initiated 
and we have found that both tensor and vector metric 
perturbations do not diverge. However a more detailed 
perturbation analysis is necessary and is complicated by 
the fact that in the presence of spatial components of a 
vector field, the decomposition theorem is violated. In 
particular, metric vector perturbations show up in the 
scalar perturbation equations [8]. However, this coupling 
may lead to observable non-Gaussianity effects in the 
CMB, and we leave this analysis for a future work '23'] . In 
closing, we believe that this model of inflation may play 
a role in connecting cosmological observables to particle 
physics in new ways, such as the baryon asymmetry in- 
dex with the measure of anisotropy and non-Gaussianity 
in the CMB. 



Appendix A: Consistency of the Analysis 

We address now the evolution and issue of backreaction 
on the 6 field due to the Chcrn-Simons term. Eventually, 
the Chcrn-Simons term could provide a back-reaction 
potential for 6 that would be responsible for the non- 
constancy of 6/Mif. However, if we use the values of 
the vector field necessary for 60-efolds of inflation we 
obtain that the maximum value of the backreaction for 
the Chern-Simons term is at the beginning of inflation. 
And in fact F^pF°'f^\^, = E ■ B\^Ja{r]o)'^ - lO'^^M^ 
is suppressed by e~^'*° at rjf. Thus the Chern-Simons 
term would provide a negligible back-reaction to 6*. It 
is straightforward to check this statement by working in 
comoving coordinates and then pull-back the result in 
conformal coordinates using dO/drj — a{t)dO/dt. From 
^ the time evolution of 6* in a homogenous FLRW met- 
ric with a — aoe^'^ is determined by 



e + iHO + 2m^e = 



E-B 

Aa^lVE 



(Al) 



We study the general solution of (jAip ignoring the Chern- 
Simons term. The conformal coupling with the metric 
provides a Hubble friction term SH9 responsible for a 
dramatic decay of 9 /Ma. in both conformal and comov- 
ing time. However, the potential acts to prevent the de- 
cay in conformal time of 9 during inflation, provided that 
m = H, namely m is order GUT scale. With an appro- 
priate choice of the initial condition, the general solution 
comoving coordinates is 6'= —6*0 e^'^*. In conformal co- 
ordinates, such a solution would give the constant value 
d9/drj — 9qH. As we have set up i? ~ M*, the require- 
ment in conformal coordinates 9/M* ~ lO^^Mp is satis- 
fied provided that 6*0 — 10~^Mp. Now we solve for 9 with 
the Chern-Simons term present, which slowly varies dur- 
ing inflation. Therefore, the particular solution to (jAip 
becomes 



_ jE-B] 
" 2almM. 



g-3//t ^ 



which at initial time of inflation contributes to the value 
of 9/M^ with a term 9p/M^ = 10~2°Mp, after washed 
out by a factor a'^{r]). And the contribution to the mini- 
mum of the field from 9p will be 9p{0)/M^ = 10~^^, thus 
negligible within respect to the general solution to (jAip . 
Therefore at initial time we find the ratio 



V{9o) ^ m^9^ 
{A-J)o {A-J)o'- 



:I0 



-10 



(A3) 



that is exponentially suppressed at later times. It follows 
that potential scalar field 9 does not directly contribute 
to driving inflation. 

Finally, we consider an analysis of the fermionic dy- 
namics in ([Ij. In our model, we argue that due to quan- 
tum squeezing, our fermions exist as a coherent state. 
Therefore, we must evaluate VEV of current operator in 
the Coherent state, which we will discover leads to the 
- redshift of the fermion current. The use of coherent 

a 

states rather than free-states is crucial in our mechanism. 
Free-states can be relevant when evaluating transition 
amplitudes between free-asymptotic particle states. In 
turn, coherent states are applied to evaluate the expec- 
tation value of quantum fields, e.g. the electric or the 
magnetic fields on flat and curved space-times, or other 
bilinear combination of fields. In our case the number 
operator is made out of with fermionic bilinears, which 
can be associated to physical observable after smearing 
out coordinate-dependence on space-time regions. Ac- 
tually, because of the time translation invariance of de- 
Sitter, smearing is performed only on space-coordinates 
and with constant test-functions. This smearing proce- 
dure must be considered as well as the expectation value 
procedure when evaluating on coherent states the classi- 
cal limit of conserved charge, such as the number opera- 
tor. 

The fermionic action has the four-fermion interaction- 
term induced by solving the torsionful part of the con- 
nection in terms of the chiral current: 



-g^^p = 



(A4) 



— 3 — — 

-itpfi) + -^(V'7'^75'0)^ + 9-0 7^6^75-0 A^ 



in which we recall that V^^ is the covariant derivative 
with respect to the Levi-Civita compatible connection 
and that y = V^e^7^. In what follows, we will use 
co-moving coordinates and recast the Dirac variation in 
terms of densitized fermions: 







-gip, V 



(A5) 



(A2) 



The equation of motion then decomposes for right- 
handed ipR and left-handed -0^ components in [4j] 

A 6 - — 

- A 6 ~ — 

«V^0L + -J jp^L^l -0 7/75-0 + q i^hAp, = ,(A6) 



8 



in which Aa = [(5* (a — 1) + 1]. In fact, the initial con- 
dition on the fermionic fields originates from the axial 
current, which is proportional to the square of fermionic 
fields. We find from the initial conditions of the currents, 
Jg^~jQ"~10^^^A/p, that the four-fermion term are negli- 
gible and easily find the solutions of (|A6p . We focus only 
on the right handed components ipR, as the solution for 
TpL will follow in a similar way. With the choice of comov- 
ing coordinates and in the Coulomb gauge, we first write 
the time component of the differential equations, namely 
^oV'k/l = c^oV'fl/L = 0. Thus, the time dependent part 
of the fermionic fields is easily recovered to be constant 
in ip, which means ip ~ 1/a^. Since fermionic fields are 
scalar functions under diffeomorphism, densitized fields 
continue to be constant in time even if recast in conformal 
coordinates. The spatial dependence of spinors, which 
turns out to contribute only for a multiplicative phase 
factor, can be found at initial time within the approx- 
imation of keeping constant the electromagnetic vector 
field. For the right-handed component we thus find 



tpR 



a 



,-»a;^ 



'-+ -I- 7 



;ioE^iAi 



e,e 



(A7) 



in which ^q and ^q are integration constants, and we have 
introduced x^ — \{xi ± 1x2). In a similar way, we find 



V'z 



xje'^o 



X§e ^^0 ^- 



AAn 



(A8) 



Notice also that the contribution from the fermionic dy- 
namics to the energy-momentum tensor, which in co- 
moving coordinates is proportional to ■0^75 Vo'0, will van- 
ish already at the classical level, consistently with our 
initial approximation Cf^={). 

Dirac fields are fundamental quantum fields and un- 
dergo second quantization through the canonical anti- 
commutation relations 

{^{x,t),ip{y,t)} = iS{x,y) 

that are imposed on space-like surfaces. As a conse- 
quence, conservation of the i/jj^tp current implies the 
constancy in time of the number operator. Indeed the 
time behavior found for the Dirac fields together with the 
plane- waves expansion of fields (on the space of the clas- 
sical solution of the Dirac equation) and the covariant 
integration procedure ensures the constancy in cosmo- 
logical time of the number operator. In order to derive 
a semiclassical quantity that could be properly associ- 
ated with the classical concept of energy-density entering 
the Einstein equations, an averaging procedure on semi- 
classical coherent states of the fermionic Fock space has 
to be considered. 

We consider a quantum mechanical system of one- 
component spinor of different flavors a, which we still 
label as V° ^-nd subject to equal-time cononical anti- 
commutation relations 

{V^Vl} = '5^ {V;",^n = WLW} = 0. (A9) 



Thus we can construct an algebra of creation ipl and an- 
nihilation ^a operators acting on |0), the vacuum state of 
the Fock fermionic space that is such that V'°|0) =0 and 
= (0|V'l ■ We further introduce Grassmann variables 
77°, which label the eigenvalues of the Ladder operators 
and whose square absolute value fjai]'^ has the meaning 
of expectation value of the number operator[4J|. As a 
consequence, we find that fermionic coherent states are 
defined by 



\v) 



— g^lv 



|0) 



— /nU'/'-''" 



and hence fulfill the relations 



r\v)=v''\v) 



\va = ivli^l 



(AlO) 



(All) 



having inner-product {r]\ri) — e''"'' . Thus, without 
a further normalization, we would find (77" 1 7/;^ 7/; 1 77) = 

This construction can be applied to the second quanti- 
zation of the Dirac fields 7/'(af, t) and ipi^j t)- Let us con- 
sider the massless Dirac field, but for facility of treatment 
let us write the Dirac equation in conformal coordinates 



{ip~e4){a^vMx,r,))=0, 



(A12) 



implying the expansion in terms of conserved space- 
momenta p and related frequency cu^ — \p\ 



a^{r))'ip{x,r]) 



E 



1 



a'^{v)^{x,v) = E 



1 



[crip)Urip)e-'P-^ + 



r p 



-yW^ 



[dip)VriP)e~^'"' 



+ cUp)Ur{p)e^P--] ^a\Tj)i;Hx,vho,{^^^) 

in which t^ is a coordinate-space normalization volume, 
the notation Urip) — ul{p)'^o and Vr{p) = vJ(p)7o has 
been meant for the spinor Ur{'p) and Vr{p), and finally 
the following orthonormality relations have been chosen 

uUp) Us{p} ^ vl{p) Vs{p) = UJp6rs , 

uUp) vsi-p) = . (A14) 



We can now proceed and quantize as in (jAlip the Fourier 
coefficient Uj^ and a- entering the expansion of fields, and 
hence define for each momentum p of the Fourier expan- 
sion a copy of the fermionic Fock space, with a vacuum 
state such that ap\Q) = for all p. Coherent states are 
therefore defined as in (jAlOp . namely 



ICp-) = e4«.^|0) , (fpl = (0|e 



p^p 



(A15) 



in which we have now hidden the flavor index a. We can 
also construct normalized coherent states |S) consider- 
ing the product |S) = Tin l^p) over each mode entering 



the Fourier expansion of the Dirac fields, and then take 
the expectation value of the number operator over such 
states. Thus the expectation value of the fermion num- 
ber density operator evaluated within the coherent state 
is given by 



n{t) = - = (S 



N 
V 






(A16) 



{E\Jd^xa^{t)^ix,t)-fo^{x,t)\E) _ no 



(S|/d3xa3(i)|S) 



m' 



in which we have introduced the total difference num- 
ber of fermions N, its density number n, and the time- 
independent part of the density number no. 

Notice that from (jA13[) it follows that Ladder oper- 
ators too must be densitized. Indeed Cr{p) ,dr{p) — > 
Cr{p}a~'^ , dr{'p)a~'^ and the same happens to their 
Hermitian conjugated operators, i.e. cl{p) ,dl{p) — >■ 
cl{p)a^^ , dl{p)a~'^. As a consequence, eigenvalues of 
inflationary coherent states are now shifted to ^(t)p = 
^pa^{t). We emphasize that inflationary coherent states 
have been defined in [30, [mI as excited states of the 
Bunch-Davies vacuum |0)bd- Thus, coherent state de- 
fined by |Sa) = ripl^p'^^) ensures that the expectation 
value of the bilinear current of the Dirac fields scales as 
a constant. Recalling indeed that J^j, = ^2(^7^7^ ejV); 
f2 summarizing the expectation value procedure on the 
state |2a), and that in conformal coordinates the in- 
verse vierbein is e^ = a^^6j, we find that J7i — 
a^"'^r2('0o757i'0o) = Ji/a, with Ji constant in time. 



Appendix B: Stability Analysis of Vector 
Perturbations 



In this model a constant Horizon size mode of the vec- 
tor fields participate initiating inflation. However, it is 
important to check whether during inflation spatial per- 
turbations in the vector fields, SAi, lead to instabilities. 
Perturbations to the spatially flat FLRW metric in con- 
formal coordinates are written as 



ds^ = a^{r]){-{l + 2'S>)dri^ -2BMx' - 
+ [{l-2'^)6ij + 2h,j]dx'dx^}, 



(Bl) 



with $ and ^ scalar, Bi vector and hij tensor Bardeen 
potentials accounting for perturbations to the metric 
tensor. It is well known that vector flelds will source 
anisotropic stress in the energy momentum tensor and 
influence the gauge invariant vector mode Bardeen vari- 
able, Bi. In this appendix we derive the equations and 
show that eventhough the vector modes in Ai are gener- 
ated during inflation, there are not instabilities because 
the solutions to perturbed Einstein equations yield a red- 
shifting vector mode in our model, _B ~ t. 



We start from 

N 



pa 



K=l 



-5]9'^'9'"d^^A%d^X^~Af^Jl^ 



K=l 

N 

and varying in A^ we obtain 

N 
K=l 



5Af^j;^^5]5A^Ji 



(B3) 



\K _ JaK 



We can further rewrite STj expressing SAi = SA^ 
diX^ (we are here following Q) 



ST;^J2^S]\SA^ J^ + {SAl<+da^)ji]-SA^j;^ + 

K=l 

- i<5;.l (5["M^9[„A^''l + d^-Af^d^JA^^^) } = 

= f:{^s)[6A^J," + {SAf + da'')J^] + 

K=l 

- i<5J^(9["M^a[„A^^'^l+9["A^a[„M^^'^l)}. (B4) 
Again, using SAf = SAf + diX^ , we find 



K=l 

+ ^ [ (^o^( ■ - %<) {d,)SA^ d,SAf^) + 
+ [dsA'< - %Af ) a[,) (Mff + 9,]x^') + 

+ (ao (%x^ + Mg) - %<5<) (a,)< - ao<) + 

- ^5] [2 {do (Mf + dsx"") - ds6A^) d[oAf^ + 
-l-9[.(a;]x^ + Mf)a[,Af]}. (B5) 

We now assume J'^ — and J'^ = -J^, with J* a con- 
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staiit three-vector, and rewrite 5T^ as 

N 



2 {do (Mf + a.x^) - dJA^) d[oAf^ + 



2^'^ 



dis {di]x'' + 5A^) d[,A^] } 



(B6) 



We now fix the gauge and impose A^ = 0, finding 



1 1 

K=l 

+ ^ [doAf^ (9.)M^ - doSAf^) + 

+ [d.Af^ - a(,Af ) a[,) (Mff + 9,]x^) + 

+ a[,(5;]X^+^A^)5[,A^]}. (B7) 

Following again [8[ , we write 



in which 



and 



n; == ^'^iVj + ('')n} + (*)n} 



'm 



d'd 



1 



5! a' a n 



'. - 3^^,- ^^ 



Thus 
SpA 



with n scalar , 

with n* transverse vector , 

(*)m 

traceless and tranverse tensor . 



ij [2 {do (Mf + a.x"") - dJA^) doA 



(B8) 
(B9) 

(BIO) 
(Bll) 
(B12) 



K=l 



d[s [d, 



h]X 



'< ' Mf d[Aif 



}+i(p^n;)<5i(Bi3) 



with 



N 



K=l 



^[doAf^ [d,)6A^ - doSAf^ 
dsAf, - a,Af ") %) (Mff + a.jx"^' 



Hj ~ "0" 



9o (af.x'" + m;^) - 9(,M^) ao4" + 



and therefore 

2 



} 



p-4n' = i:{-^(aa"+X4f)j-+ 

+ Jj [aoAf (aM^ ~ SoMf ) + 



(B14) 



(B15) 



+ 



(ao (ax"" + -^^f ) - 5,^0 aoAf + 



+ a^, ( 6Af + 9,]x^) a[,Aff ] } 



^*] 



TV 



E{-^(^»x''+^^f)-^^ 



4 



a 



■ [9oAf (a<5<-ax''-'Jif) + 

+ a[,Mfa[,Aff]}. 

'"'Dj components of the perturbed 
energy-momentum tensor 

We want here to obtain the '^"'ffi components of the 
perturbed energy-momentum tensor H' and prove that 
(^)n*- = 0. We consider then 



N f 
A'=l '- 



K\ ji 



X'^ + 5Af J 



doAf (d,5A^ - ax^ - 5Af) + a[,M^ a[,Aff 



} 



and focus on its components involving di derivatives. For 
simplicity, we work at K fixed 



ao^f (d,5A^ - d.x" - 5Af) + d^s5Af^ d^.A 
-f(ax^ + Mf)J* 



(B16) 



and analyze each term involved in this latter equation. 



• About the first term in (|B16p we notice that, be- 
cause of the constancy of J* , we can rewrite it as 



11 



di (x^J')- But if this term belong to ^'"^U^, then 
we should have 



d.{x''r)=o.. 



(B17) 



which happens iff diX^ = 0, i.e. if the term is not 
present at all and 6Af does not include perturba- 
tions coming from the gradient of a set of scalar 
functions x^ ■ We can conclude that this term can 
not contribute to ("^11'. 

• We apply the same reasoning to the third and the 
forth term, assuming the transverse Coulomb gauge 
for A^, namely d, Af = 0. Then 

if a. (Mf - xn = 



d, 
d^ 



AfiSA^^x"") - (9.Af ) (Mf - x^) = 



Af {6A^ - x^) 



(B18) 



If this term is the contraction of the i and j indices 
of term belonging to '^"•'ll*, then it has to be zero 

by definition of *^^)n*, i.e. 



Af(Mf-x'^) =0 iff SA^^x""- (B19) 



But this condition would imply that the term which 
should contribute to ^"^II* is identically zero. 

• The remaining term involving derivative in di that 
we should consider is the last one of (IB16|) . We can 
rewrite this last term, again using di Af = 0, as 



d[s5Af^d^iA 



\K 






(a[,Af ) Mf 



-dAduAfASA"^ = 



-a,(a,Af-a,Af)Mf.(B20) 



We focus now on the last term in (|B20I) . Using 
again di Af = 0, we find 



-a, {dsAf - a,Af ) = -a, aaf 



^^^^A^ = 



(B21) 



where we have used in the last hand side the equa- 
tion of motion for A^ . Thus this last contribution 



in (JB20I) can not result to be the trace in i and j 
of a term belonging to ^"^11*. 

The analysis of the points above ensure us that '^"•'ll* = 0. 
As a consequence, within the same notation of [8| , we find 
for the vector perturbations B^ that 

djB' + 2- djB' = 87rGa2p^(")m = , (B22) 

a ■' 

implying 






(B23) 
The differential equation regulating the evolution of the 
transverse vector perturbations 5A^ reads 

dl5A, - djd' 5A, ~ doA^d^B, ~ 2doh,jdnA^ - 5* = 

(B24) 
in which Si is defined by means of the decomposition of 

2<i>d^A^ -I- ao($ + *)aoA; = ^^S + S, (B25) 

in a scalar S and Si, which is a transverse vector. Notice 
that, assuming $ = -^, (|B25p implies S ~ e^'^'^ f [t]) , 
with f{rj) generic. When $ = 0, then both 5i and 5 
are vanishing, and the only vector perturbations Bi and 
tensor perturbations /ly enter (jB24|l . These latter are 
determined by 



daa 



Wi 



ao^/iy +2-^ao/i.y -a*afe/i,j = s^Gp^'^ny ~ o , (b26) 

in which we are disregarding the anisotropic stress ^^'liij. 
The solution to ^^ is given by h,, ~ ^----m^/<n)) ^ 
Thus, the evolution of (jB24p is sourced by terms that are 
suppressed in a{ri), and the stability of the solution for 
Ai{7]), and therefore a{rj), is recovered. 
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We thank Jim Gates for pointing out this fact to us. 
It is possible that this scalar field is not fundamental and 
is instead analogous to a gauge singlet quark condensate 
such as the rj meson. It is well known that the rj particle 
couples to the Abelian Chern-Simons term to address the 
(7(1) problem. 

We will be justified to not include the coupling of the 
superhorizon to the vector modes Vi, since we will show 
that spatial anisotropies on superhorizon scales will be 
suppressed by \/iV. We thank Justin Khoury for raising 
this point. 

Even if Pk (k > 9/M,) is imaginary and the expo- 
nential of the general solution in @ and ^ become 
oscillatory, the gauge field still grows because of the 
J-/[2Il\l-IIr])\ terminH[i/,/3, J_,77], due to the back- 
reaction provided by the gravitational field. 
It has been shown by Ford 2| that vector fields neces- 
sarily generate an anisotropic metric arising from off- 
diagonal components of the energy momentum tensor. 
The anisotropic parts of the energy-momentum tensor 
due to the field-strenght are Tot = Si, where S = E x B 
is the Poynting vector, and Tij = EtEj -\- BiBj — [E^ -\- 
B'^)5ij/2. Likewise, we will find that values of the gauge 
field necessary to generate infiation in the isotropic part 
of the energy-momentum tensor yield electric and mag- 
netic components that are perturbations. This occurs be- 
cause the derivatives that act on A to generate the elec- 
tric and magnetic fields are suppressed by a~*. Due to 
the growth of the gauge field the electric field will redshift 
as a~^. We thank Robert Brandenberger for bringing this 
to our attention. 

Numerical simulations show that A- J is monotonically 
decreasing in time during infiation and has a maximum 
relative variation of ~ 0.5. 

We choose the initial conditions A*^ = A^ — 0. 
For more details on the role of the Chern-Simons term in 
cosmology we refer the reader to [19l - [2ll ]. 
We thank the referee for raising this crucial possibility. 
We are grateful to the referee for raising the point of 
non-standard charge assigments in our model. 
Since we have extended the gauge sector to include A'' 
copies of gauge field, we are assuming that these abelian 
factors only carry solely leptonic charge in the minimal 
standard model. In general we have the freedom to couple 
the gauge fields to other states that are not charged under 
the weak-hypercharge assignment (such as a dark sector), 
but this falls outside the scope of this present work. We 
thank the referee for raising this crucial point, since non- 
standard assignments may provide a mechanism of dark- 
matter production at the end of infiation. 
In general, fermionic coupling to general relativity in- 
duces torsion. Inclusion of a four-fermion term will elimi- 
nate the torsion and yield metric compatibilty. However, 
the contribution to the four-fermion interaction is negli- 
gible, as it is proportional to J^ , whose components are 



13 

suppressed by a factor of a. Furthermore, the introduc- namely \0)ri°' = r^^lO) and \0)fj'^ — 77"|0). The same prop- 

tion of densitized Dirac fields will introduce an additional erty holds for {0|. Moreover, we can say that jy" is the 

redshift factor proportional to y^^-g. complex conjugate Grassmann variable associated to the 

[44] Notice that, being 77" a Grassmann variable, we have that Grassmann variable r?". 
the vacuum states commutes with Grassmann variables. 



